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ABSTRACT 

A  flat  strip  antenna  Is  embedded  In  a  planar  dielectric  slab,  which  Is 
surrounded  on  both  sides  by  layers  of  compressible  Isotropic  electron  plasma. 
Several  closed  form  expressions  are  obtained  for  the  impedance  with  a  sinusoidal 
current  distribution  along  the  antenna.  Die  antenna  impedance  Is  computed 
numerically  when  considering  the  perturbations  of  the  antenna  current  by  surface 
waves.  Except  for  thin  plasma  layers,  the  antenna  impedance  can  be  computed 
using  the  same  current  distribution  as  for  a  cold  plasma.  This  supports  the 
validity  of  earlier  work  which  neglects  the  perturbation  of  the  antenna  current 
by  plasma  waves.  However,  It  Is  essential  to  consider  the  finite  transverse 
dimension  of  the  antenna  and  the  presence  of  an  Insulating  layer  or  of  an  Ion 
sheath. 
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The  variational  calculation  of  the  admittance  of  a  slot  antenna  Immersed  In 
a  two-layer  compressible  plasma  was  recently  reported  by  Galejs  [1966a]  and  this 
paper  contains  also  a  detailed  discussion  of  past  work  in  tills  problem  area.  The 
surface  waves  supported  by  the  plasma  layers  were  shown  to  modify  the  field 
distribution  along  the  antenna  and  to  affect  the  antenna  admittance.  A  similar 
effect  can  be  anticipated  also  for  a  finite  linear  antenna,  which  Is  treated  in 
the  present  paper. 

A  flat  strip  antenna,  shown  In  Fig.  1,  Is  embedded  in  a  planar  dielectric 

slab  of  thickness  2A,  which  is  surrounded  on  both  sides  by  a  compressible 

isotropic  plasma  of  thickness  h.  The  dielectric  layer  is  Intended  to  approximate 

the  effects  of  the  Ion  sheath  which  Is  formed  around  the  antenna.  Die  expressions 

for  the  fieldsare  derived  starting  out  with  linearized  field  equations  [Oster,  i960] 

as  Indicated  In  Section  2.  The  antenna  impedance  is  formulated  in  Section  3  In 

the  same  way  as  for  a  linear  antenna  In  an  Incompressible  plasma  or  dielectric 

layers  [Galejs,  1965b].  The  antenna  current  distribution  is  assumed  to  be  the  sum 

of  two  terms.  The  first  term  is  a  sine  wave  with  the  same  wave  number  k.  as  In 

A 

incompressible  plasma.  The  second  term  is  representative  of  surface  waves  sup¬ 
ported  by  this  plasma  geometry.  The  wave  number  kg  of  the  surface  wave  Is 
determined  from  the  solution  of  the  transcendental  equation,  which  describes  the 
poles  of  the  Integrands  of  the  field  expressions.  The  resulting  Impedance  express¬ 
ions  are  evaluated  numerically,  as  is  done  for  an  incompressible  plasma,  except 
that  the  integrals  consider  In  detail  the  contribution  of  the  surface  wave  poles. 

In  the  special  case  of  an  antenna  Immersed  in  a  homogeneous  compressible  plasma 
with  an  assumed  sinusoidal  current  distribution,  closed  form  approximations  of 
the  antenna  impedance  are  derived  In  Section  5,  which  compare  closely  with  the 
results  of  other  investigations  [Seshadri,  1965;  Balmain,  1965 ] .  Burner i cal 
results  are  discussed  in  Section  6. 


2.  FIELD  EXPRESSIONS 

•Ihe  field  equations  for  a  suppressed  exp(-iwt)  tine  ‘variation  can  be  written 


following  Oster  [i960]  as 

7  x  E  a  iuji  H 
—  0  — 

(1) 

7  x  H  «  -  iw€Q  E  +  N  e  V 

(2) 

(  -  i«  +  v)  m  HV  ■  +  NeE  -  7  p 

(3) 

u2m  R  7.V  -  iwp 

(k) 

where  E  and  H  are  the  electric  and  Magnetic  field  vectors,  I,  V,  v  are  the 
average  values  of  particle  density,  velocity  and  collision  frequency  respectively; 
p,  u,  e,  and  m  are  the  tine  varying  component  of  the  scalar  pressure,  the  acoustic 
velocity  and  particle  charge  and  mass  respectively. 

Eliminating  V  from  (2)  and  (3)  gives 


E  »  - 

------  7  X  H  -  7  p 

iwt  «  -  tie  * 

0 

(5) 

with 

c  -  1 

w(w+iv) 

(6) 

2 

(1)  at 

p 

n  e2/(«Qm) 

(7) 

Substituting  (5)  in  (l)  results  in 

^H+k2H-0 
**  0  ** 

(8) 

2  2 

with  k  ■  w  u  f  «.  For  the  later  work  it  is  convenient  to  separate  the  vector 
0  0  0 

components  of  H  into  the  TE  and  B!  parts  that  can  be  derived  from  two  scalar 
functions  ¥  and  •  each  of  which  satisfies  a  wave  equation  similar  to  (8).  H  , 
H  and  H  are  related  to  ¥  and  •  by 

y  * 


H*  ’  ikr  (  m *  +  4  *  ) 

(9) 

Hy  *  l5T  (  Say  1  -  V2  If  *  ) 

(10) 

(11) 

2 


Eliminating  E  from  (2)  and  (3)  gives 


1--  sH  [ 


V  x  H  +  ii  V  p 


(12) 


and  substituting  (12)  in  (4)  it  follows  that 

+  k  2  p  =  0 
P 

with  k  2  »  [w(u*iv)  -  w  2]/u2  . 

Xr  Ir 


The  solutions  Y,  #  and  p  for  the  three  regions  of  Fig.  1  can  be  expressed  in 
terns  of  the  integrals 


f  f  1  (v).-*- ♦  l.V) 


du  dv 


-to 


to  » 

f  f  e‘lTy  (.r“‘  +  V"7™’)  4u  dT 

CO  -to 
CO  CO 

f  f  Cn(*,v)e-im  .-l7y  ♦  8cn  e‘V)  du  dv 


where  the  subscript  n  «  1  or  2  designates  the  region,  the  amplitude  *  0, 


the  reflection  coefficients  R  _  =  R.  -  *  0  and  where  7  «  i  vk  *  -  u  -  v  and 

a  j  d  j  en  en 

x  g  2  2 

7pn  *  i  vkpn  -  u  -  v  are  defined  in  the  second  quadrant  of  the  complex  plane 

including  the  positive  Imaginary  and  the  negative  real  axes .  Also,  k  -  k  „  ■ 

ei  ej 

k©  =  « •J uq«o  .  After  computing  the  components  of  H  from  (l^)  and  (15)  with  the 
aid  of  (9)  to  (ll),  the  components  of  E  and  V  are  obtained  from  (5)  and  (12)  by 
making  use  of  (l6). 

In  previous  work  [Galejs,  1965b]  a  stationary  expression  for  the  antenna 
impedance  was  derived  in  terms  of  the  electric  fields  and  currents  of  the  antenna 
which  were  related  to  the  coefficients  A^  and  of  the  functions  Y^  and  for 
the  region  just  outside  the  aperture.  In  the  present  problem  of  a  compressible 
plasma  the  fields  can  be  specified  with  the  sne  coefficients  A^  and  after 
relating  them  to  the  remaining  constants  R^,  R^,  *2.  V  V  V  V  V  *3 


3 


and  B..  of  the  field  expressions  (14)  to  (l6)  from  a  total  of  IX)  boundary  conditions 
(E,E,H,H,  continuous  and  V  ■  0  at  2  «  t.  and  t.).  These  equations  can  be 
solved  for  the  reflection  coefficients  and  .  A  lengthy  algebraic  nani- 
pulation  results  in 


\i  *  ' 


/  Z\2ll 

27elzl  7.x(8  +  Se2 


^  V  ^  -  Jk  1 

p  ♦  * )  * 7<2  ^  i ) 


’’ell8 


R»2  *  8 


2re2*2  7«2  '  7e3 
7.2  +  7.3 


with 


^1*1  P  -  Q 
8  P  +  Q 


<-£  ( •** - v)  -  v  ( ♦ * ) 


27^*1 


“7e2h  Z'  27e2Z2 
+  d^Cs  e  ed  (  e  Bt£  d  + 


**) 


(21) 


k-2  \2  ‘7e2h  7e2 

-  )  +  4,Ct  -  d,Cs  e  8<:  +  — 

3  )  3  3  7e3 

—  'j  -  d,Ct  +  4,  Ob  1*^  + 

-o  )  3  3  7.o 


27e2®2 


•(5J- 

(M- 


where 


(1  -  ^)(»2  +  v2) 

7.J  V 


Ct  -  coth  (7p2^)  (2U) 

Cs  =  l/sinh(7p2h)  (25) 

The  reflection  coefficient  of  the  TE  nodes  R  ,  of  (17)  and  (l6)  are  the  sane  as 
for  an  incompressible  plasma.  [Galejs,  1965b].  The  compressibility  of  the  plasna 


l 


introduces  the  terms  proportional  to  dj  in  the  expressions  (l$)  and  (22)  focr  the 
reflection  coefficients  R^j  of  the  TM  modes.  R^  of  (19)  is  the  reflection  co¬ 
efficient  for  radiation  in  a  compressible  plasma  layer  of  thickness  h.  In  the 
limit  of  h  -*  co,  Rfel  »  0  and  Rfa2  of  (22)  is  the  reflection  coefficient  for  radiation 
into  a  plasma  halfspace  vhich  is  separated  from  the  antenna  by  an  ion  sheath  of 
thickness  ■  A.  These  two  limiting  forms  of  reflection  coefficients  of  the  TM 
modes  have  been  obtained  previously  [Galejs,  1966a].  It  can  be  also  shown  that 
R&1  and  R^  0  as  h  -*  0,  if  kgl  =  k^. 

The  assumption  of  a  rigid  boundary  of  the  plasma  layers  with  V  =  0  can  be 

z 

replaced  by  more  elaborate  boundary  conditions  of  continuous  scalar  pressure  p 

and  continuous  normal  particle  velocity  V  .  This  leads  to  increased  algebraic 

z 

complexity  of  the  derivations,  but  give  essentially  the  same  results  as  the  simpler 

and  less  accurate  boundary  condition  of  V  »  0  [Galejs,  1966a] . 

z 

3.0  AHTEHNA  IMPEDANCE  AND  RADIATION  RESISTANCE 


The  ion  sheath  surrounding  the  antenna  is  approximated  by  free  space  and  the 
antenna  Impedance  can  be  formulated  in  terms  of  the  TE  and  TM  modes  of  this  region 
as  for  a  linear  antenna  in  a  stratified  dielectric  [Galejs,  1965b].  The  surface 
current  of  the  antenna  is  assumed  to  have  only  an  x-component  J^,  vhich  in 
addition,  is  an  even  function  about  x«0,  y«0  and  the  driving  point  Impedance  of 


a  flat  6trip  antenna  may  be  computed  from  the  expression 


[2*  l(x*0)r 


//  F(u'T)  [  //Jx(x-y) 


C08  UX  COS  vy  dx  dy^  (26) 


where 

F(u,t) 


(l/a)  r  vf_  _  “2  7 el  1~Bbl 

»a+va  L7el  '  kel2  I4«bl 


(27) 


The  antenna  impedance  Z  depends  on  the  current  distribution  Jx(x,y)  along  the 


antenna,  vhich  is  assumed  to  be  in  the  form 
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Jx(x,y)  *  |  A  3in[kA(j-|x|)]  +  B  sinCkgCj-lxj)]  j-  f(y)  (28) 

with  f(y)  =  l/(2 «)  »  const.  The  phase  constant  kA  is  selected  in  such  a  way  that 
the  impedance  Z  computed  with  B=0  approximates  the  impedance  Z  of  the  correspond¬ 
ing  cold  plasma  geometry,  which  was  derived  in  an  earlier  formulation  using  a 
more  accurate  two  term  trial  function  [Galejs,  1965b].  The  second  term  of  the 
current  distribution  (28)  should  approximate  the  effects  of  the  finite  plasma 
temperature.  The  phase  constant  k^  is  representative  of  surface  waves  that  are 
supported  by  the  antenna  in  the  presence  of  the  compressible  plasma  and  it  will 
be  discussed  in  Section  4. 

Substituting  (28)  in  (26)  and  utilizing  the  stationary  character  of  the 
Impedance  expression  for  determining  the  complex  amplitudes  A  and  B,  it  follows 
that  [Galejs,  1965b] 


Z  * 

with 

A  = 


7BB  FA  ~  yAB  FB 
7aa  fb  "  7AB  7 A 

7AA  7BB  “  7AB 


*B  7AA  ^A  FB  7AB  +  *A  7HB 


(29) 

(30) 

(31) 


(32) 


4iwp  r  /  .  \2  r 

7NM  *  -J*  V)  / 

o  o 

h 

gjj(u)  =  -A- "2  (cos  Ui  -  COS  kjjl) 

Vu 


(33) 

(34) 


Ihe  impedance  can  be  computed  also  for  the  sinusoidal  field  distribution 


term  by  setting  B 


0  in  (28).  This  gives 


(35) 
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4.0  SURFACE  WAVES 


The  phase  constant  k_  is  assumed  to  be  the  same  as  for  surface  waves  which 

rs — 9 

are  supported  by  the  plasma  layers,  and  k^  »  vu'  +  v*  Is  determined  from  the 
poles  of  the  integrand  (26),  where 


(36) 


Substituting  (19)  to  (22)  in  (36)  results  in 
/  K o  \2  7.„  /  k.o  \2 


WJkY  !k  ,  .  /.  .  %a\/ 

\l  '  W  ‘  re3  '  kel  )  el  1  ^  \3  3  Vl 


coehre2h 


-  Ct 


+  tanh  re2h  {  [(  ^  )  -  a3«]  [(  ^  )  coth  rel«x  +  ^ot] 

-  jf  if  +  dld3(Cs)2  }  -  0  (37) 

where  the  symbols  Cs  and  Ct  are  defined  in  (24)  and  (25).  In  the  limiting  case 
of  h  -*■  0  and  kj  «  kp2,  (38)  is  reduced  to 

ootbr^i  -  ^  (38) 

Equation  (38)  can  be  seen  to  be  identical  with  (26)  of  Oalejs  [1965b],  which 
denotes  the  characteristic  values  of  urface  waves  for  TM  modes  in  a  dielectric 
layer  of  relative  dielectric  constant  *r  »  (k^/k^)  and  of  total  thickness  2Zy 
Introducing  infinitesimal  plasma  losses  (v/o)  in  the  limit  of  h-*«  , 
tanh  7egh  and  Ct  -*•  (-1)  and  Cs  -*>  0.  Equation  (37)  simplifies  now  to 


coth  y  ,  z.  -  d,  +  - 

el  1  1  y  , 

el 


(39) 


which  can  be  also  obtained  from  (36)  and  (19)  by  setting  R^  *  0  in  (20)  and  (2l) 
Equation  (39)  can  be  expected  to  have  solutions  for  k^  *»  k^2  »  k^j  .  Letting 


«2  =  (kg2/kel)  (39)  can  be  expressed  as 


7 


1  -  €, 


«2  coth  k^  = 


■J 1  -  vv 


For  thick  insulating  layers,  kgZ1  »  1  (40)  is  changed  to 


B  p2  ~f=  ^ 

2 

which  represents  surface  waves  guided  along  the  boundary  of  a  seal -Infinite  plasma 
(Eq.  (54),  Hessel  et  al.  [1962];  Eq.  (55) »  Oalejs  [1966a]).  For  thin  insulating 
layers,  kg  «  1  (40)  becomes 


*b-v[1  +  5  ( Vi)2  +  •••  ] 


and  the  wave  number  of  the  surface  wave  k_  approaches  the  wave  number  k  .of  plasma 

o  J6 

waves  in  a  homogeneous  plasma.  In  the  limit  of  «2  -*■  0  and  tg  «  kfiz1  both  ( 1+1 ) 
and  ( U2)  are  replaced  by 


kB-ko  k 


where  k  *uv(i  (  .  A  more  detailed  examination  of  ( 39)  shows  that  ( 40 ) 

000 

represents  its  only  solution  and  that  kg  lies  within  the  limits  indicated  by  (4-1) 
to  (43)  if  <2  is  in  the  range  0  <  *2  <  1. 

For  a  plasma  layer  of  finite  thickness  ’h’  (37)  will  be  simplified  by  assum¬ 
ing  kg  *»  kp2  »  kgj,  and  h  of  the  order  of  the  free  space  wavelength  2«/kQ. 

This  gives 


6 


For  h  large  and  |kg|  >  |kp2|  the  right-hand  side  of  (44)  approaches  aero.  The 
zero  of  the  first  square  bracket  represents  a  wave  guided  along  the  plasna  bound¬ 
ary  near  the  antenna  ( 1*0 ),  while  the  zero  of  the  second  square  bracket  represents 
kg  of  a  wave  guided  along  the  outer  plasma  surface  which  is  the  same  as  kg  of  the 
wave  guided  along  the  inner  boundary,  if  kgZ^  »  1.  These  two  waves  are  coupled, 
if  the  right-hand  side  of  (44)  is  finite,  but  the  resulting  value  of  kg  will 
approximate  kg  computed  from  ( ho) .  Equation  (44)  has  further  solutions  with 
|  kg |  <  |  k^  | .  However,  it  can  be  shown  that  these  waves  are  excited  with  small 
amplitudes  if  the  plasma  layer  has  a  thickness  comparable  to  the  free  space 
wavelength  [Galejs,  1966b].  These  waves  are  not  considered  at  present. 

A  real  root  kg  is  determined  first  from  the  solution  of  (44)  for  a  lossless 
plasma.  For  a  lossy  plasma  the  complex  root  is  found  by  applying  Newton's 
iteration  method  and  by  using  the  former  reed  root  for  the  initial  estimate. 

For  h  »  z1  the  roots  of  (44)  differ  negligibly  from  the  roots  of  (40),  which 
lie  in  the  range  indicated  by  (4l)  to  (43). 


5.  IMPEDANCE  OF  A  SHORT  ANTEKNA  IN  AN  UNBOUNDED  PLASMA 

The  impedance  of  a  short  antenna  (kQi  «  l)  will  be  computed  from  (35)  for 
an  unbounded  plasma  with  an  assumed  sinusoidal  current  distribution  of  kA  »  kg2  . 
For  z^  =  0  and  h  -♦  «  *  Rg2  ■  0,  and  the  substitution  of  ( IT)  and  (19)  in 

(27)  shows  that 


F(u,v)  = 


2(u2  +  v2) 


( 


.2 

e2 


u27 


e2 


e2 


) 


2 


u 


l-«. 


2  k  _  7  o 

e2  p2 


(45) 


The  first  term  of  (45)  does  not  depend  on  the  finite  temperature  of  the  plasma, 
and  substituting  it  in  (33)  the  impedance  dee  to  the  electromagnetic  waves 
(subscript  e)  can  be  computed  from 


00 


Z  =  R  +  i  X  = 
e  e  e 


2  *^0 


f 


du 


(cos  uf-cos  kAf)‘ 

32  2 

k.  -  u 
A 


[»( 


8incv 

tv 


(46) 


e2 
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The  antenna  resistance  Is  computed  by  confining  the  range  of  integration  to 

2  2  2 

values  of  u  and  v  where  u  +  v  <  k.  .  Over  this  range  of  integration  the  sine 

A 

and  cosine  functions  can  be  replaced  by  their  small  angle  approximations .  The 


integrations  are  elementary  and 

He  -  ao  ^  (V)2 


The  reactance  is  computed  as 


2u“o 


09  00 


(cos  ui  -  cos  kAi)‘ 


r  p  p  r  r  i  'cos  ^  “ 

[  /du  J  dv  +  J*>Jto\  -J* 

ri  — — —  n  '  i  ' 


+v  -k , 


sin  cv 

«v 


)2  (*») 


4E? 

A 


The  v-lntegrations  become  elementary  by  noting  that  sin  cv  *=  cv  over  the  range 

r  p  2 

of  v  where  vu  +  v  ^  v  ,  which  applies  strictly  if  i  -»  0  and  u  <  u  is  finite. 


The  v  integrals  are  evaluated  to  give 


2  up 
o 

(tk/V 


(cosui-cos  kAi)' 

7T2  2\ 

(kA  -  u  ) 


!-c'log(«kA)  •  |  log  |l  -  (  i)2| 

A  w 


where  C  =  0.5772... is  Euler's  constant.  The  integral  which  is  proportional  to  the 
constant  term  of  the  square  brackets  is  evaluated  directly.  The  remaining  term 
of  the  integral  which  Involves  a  logarithmic  function  is  reduced  to  tabulated 
Integrals  after  expanding  the  denominator  in  partial  fractions  and  changing  the 
variables  of  integration  to  y  =  (ix/k.)  +  1.  This  results  in 


Xe=^7  [log  ^  +  0.5  -2  logs] 


The  contribution  of  the  plasma  waves  (subscript  p)  to  the  antenna  impedance 
is  computed  using  the  second  term  of  (U5) .  This  results  in 


R  +  iX 
P  P 


2l*o  1-*s  r 

*2/2  «„kS 1 

d  O  O 


COSU/-] 


up 

u 


slncv 


K-  'K  * 


The  computations  will  be  made  first  for  very  thin  antennas,  (k^O  «  1. 
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*  JH-V 

M 


2  2  2 

In  the  resistance  calculations  u  +  v  <  k^  and  sine  «v  •»  «r  over  this  range  of 


integration.  It  follows  that 
120(1- 


P  *2*o  1 


—  {a  sidc^/)  -  Si(2  k^t)  -  --Bl-°^fg^g)  |  (52) 


For  kpgf  »  1  (52)  simplifies  to 

V  W 

and  for  k^f  «  1  (52)  is  approximated  by 

RP  - 10  d-V  (V>2 


(53) 


(54) 


Also,  Rp  of  (54)  can  be  obtained  directly  from  (51)  by  using  snail  argument  approxi- 

2  2  2 

nations  of  the  cosine  function.  In  the  reactance  computations  u  +  v  >  k_  . 

p2 

The  v-lntegrals  are  evaluated  using  the  same  approximations  as  in  (48),  which 
gives 

2«40a-^) 


*2(’*o<)‘ 


■f du  <1-eo|  HOf.  1 . 

n  ^  — 


'll 

1  kp2  J 

(55) 


The  Integral  which  is  proportional  to  the  constant  term  of  the  square  brackets  is 
evaluated  directly.  The  remaining  ten  of  the  Integral  which  Involves  a  logarithmic 
function  is  evaluated  using  the  relation 

22, 


®  2  2  p 

f  dx  lo.g -Lgl* JL-L  (cos  xy2  -  cos  xy^  -  x  J  Cl  (|)  dy 


(56) 


where  Ci(x)  is  the  cosine  integral.  (Equation  (56)  is  obtained  by  integrating  a 
tabulated  pair  of  Fourier  sine  transforms  [Oberhettinger,  19571.)  For  k^f  »  1, 
(55)  becomes 

120(1-«J 


XP  ~  " 


UU-«oJ  pa  1  r  1  1  1 

rrr-  {f  - c  -  V  -  77^  |_2  coe(V>  -  s 606  (2V>J  }  ™ 
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For  antennae  which  are  vide  relative  to  the  plasma  wavelength,  k^c  »  1, 

2  2  2 

only  those  values  of  u  and  v  where  u  +  v  «  kpg  will  contribute  significantly 

to  the  integrals  and  u  and  v  can  be  neglected  relative  to  k^g.  The  integrations 

become  elementary  in  (51)  and 
60jt  (1-«J 

Z  *  R  =  — r— r— (58) 

P  P  *2*0*  *P2 

It  may  be  noted  that  (54)  has  been  obtained  first  by  Hessel  and  Shmoys  [1962], 
and  (52)  has  also  been  derived  by  Seshadrl  [1965]*  Balmain  [1965]  has  considered 
a  cylindrical  dipole  antenna  of  a  finite  radius  r,  and  expresses  the  antenna 
Impedance  in  terms  of  Bessel  functions.  For  r  small  his  resistance  agrees  with 
(53),  and  the  reactance  is  the  same  as  the  leading  term  of  (57)  if  the  term  (3/2) 
is  replaced  by  2  log  2.  However,  there  are  differences  for  r  large.  His  resistance 
is  proportional  to  cos  (k^r  -  n/k)  with  a  peak  value  4  times  larger  than  shown 
in  (58).  His  reactance  can  have  either  sign  and  has  a  peak  value  two  times  larger 
than  R  of  (58).  However  the  reactance  X  of  the  present  calculations  is  negligible 

Jr  X* 

relative  to  R  for  k  _e,  »  1. 

P  P2  1 


6.  DISCUSSION  OF  NUMERICAL  RESULTS 


The  antenna  Impedance  has  been  computed  for  very  thick  plasma  layers  (h  -*•  ») . 
The  antenna  has  a  finite  width  (fl  =  2  log  4i/e  -  10)  and  the  medium  has  slight 
losses  (tan  ftg  «  0.03).  The  calculations  are  made  for  very  thin  (A  -  10~^x)  and 
for  (A  »  10”2x)  layers  of  insulation  and  the  results  for  a  warm  plasma  with 
c/ug  «  100  are  compared  with  a  cold  plasma  (c/Ug  =  »).  For  small  values  of  Cg 
the  antenna  is  thin  relative  to  the  acoustic  wavelength  (kpgC  «  0.43  for  «g  =*  0.01) 
and  the  appropriate  closed  form  approximations  are  indicated  by  heavy  lines  in 
Fig.  2.  The  finite  plasma  temperature  tends  to  increase  the  antenna  resistance  and 
to  decrease  the  capacitive  reactance.  Also  the  net  change  due  to  the  finite  plasma 
temperature  appears  to  be  larger  for  small  values  of  «g. 

Landau  damping  of  the  longitudinal  plasma  oscillations  may  become  significant 
if  kfi  becomes  comparable  to  upg/ug.  Noting  that  kp2  “"kg  ,  this  damping  will 


12 


be  small  if  kp2  <  w^/u^  .  This  condition  can  be  rearranged  into  w  <•1?  .. 
or  «2  <  0.5.  The  damping  can  be  neglected  for  the  smaller  values  of  «g  shewn 
in  Pig.  2. 

The  further  numerical  calculations  are  restricted  to  a  constant  value  of 
c2  *  0.03  and  the  variation  of  the  antenna  impedance  with  the  velocity  ratio 
c/ug  is  shown  in  Pig.  3.  The  closed  form  approximations  of  k^t  «  1  and  »  1 
apply  to  c/ug  small  and  large  respectively  and  are  indicated  by  heavy  lines. 

The  finite  width  of  the  antenna  can  be  neglected  only  if  k^e  «  1.  Approximation 
of  a  finite  antenna  by  a  current  filament  will  Introduce  considerable  errors  if 
kp2€  »  1,  which  corresponds  to  c/u2  >  1000  in  the  present  numerical  examples. 

The  data  shown  in  Figs.  2  and  3  refer  to  an  antenna  of  a  constant  length 

/  *  0.25X  and  the  impedance  Z  of  ( 30)  computed  with  a  two  term  trial  function 

differs  by  less  than  a  few  percent  from  Z  of  (35)  calculated  for  the  sinusoidal 

6 

field  distribution.  Further  investigations  show  that  the  current  ripple  with 
wave  number  kg  has  practically  no  effect  on  the  antenna  Impedance  if  the  antenna 
length  2 /  is  an  integer  multiple  of  the  wavelength  of  the  surface  waves  or  if 
Re  kg/  *  nn,  where  n  is  an  integer,  The  plasma  wave  components  of  the  current 
distribution  cause  a  maximum  perturbation  of  the  antenna  impedance  if  Re  kg/  ■ 

(n  +  0.5) a.  For  a  lossless  plasma  the  wavenumber  kg  is  real,  and  V  «  nit 
corresponds  to  a  standing  wave  pattern  with  a  null  at  the  antenna  center.  Such 
a  standing  wave  does  not  perturb  the  antenna  current  at  the  feed  point  and  will 
not  be  excited.  Similarly  kg/  =  (n  +  0.5)n  causes  a  standing  wave  pattern  with 
a  maximum  at  the  antenna  center,  it  will  perturb  the  antenna  current  at  the  feed 
point  and  will  affect  the  antenna  impedance. 

The  antenna  Impedance  is  shown  in  Fig.  4  for  a  two-term  trial  function  (z) 
and  for  an  assumed  sinusoidal  current  distribution  (Z  )  for  vaiious  amounts  of 

D 

plasma  losses.  The  antenna  length  is  made  equal  to  /  -  8.5it/Re  kg.  Both  sets 
of  computations  give  nearly  the  same  resistance  for  tan  &2  <  0.1.  The  amplitude 
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ratio  A/B  of  the  tvo  components  of  the  antenna  current  distribution  Is  computed 
from  (29)*  It  is  nearly  500  in  magnitude  for  small  values  of  tan&2  and  is 
approximately  160  for  tanfig  ■  0.3. 

The  effects  of  gradually  Increasing  the  thickness  of  the  antenna  insulation 
are  shown  in  Fig.  5.  Hie  antenna  resistance  is  decreased  by  Increasing  thickness 
of  the  insulation  A  .  However  the  oscillatory  peaks  of  the  antenna  resistance 
which  are  observed  for  larger  values  of  A  are  observed  also  for  a  cold  plasma. 

For  small  values  of  A  the  antenna  resistance  in  the  cold  plasma  case  is  due 
principally  to  plasma  losses  and  R  «  tanbg  X.  The  Impedance  data  shown  in  Fig.  5 
has  been  computed  for  an  assumed  sinusoidal  current  distribution  from  (35)  and 
the  use  of  a  two-tem  trial  function  gives  a  change  of  less  than  5  percent  in  the 
impedance  figures. 

Antenna  impedance  for  plasma  layers  of  various  thicknesses  h  is  shown  in 
Fig.  6.  The  impedances  Z=R  +  iX  and  Z  *  R  +  iX  are  both  indicated.  For  h  >  0.3X 

the  antenna  impedance  is  nearly  the  same  as  for  h  =  °°  .  However  for  small 

values  of  h  the  antenna  impedance  exhibits  several  high  resonance  peaks.  The 

largest  peak  occurs  near  h  =  0.028X  where  k  _h  «  it.  These  peaks  may  be  attributed 

P<= 

to  the  resonance  of  the  plasma  waves  in  the  transvers  dimension  of  the  plasma 
slab.  The  surface  waves  which  are  guided  along  the  plasma  slab  have  a  consider¬ 
able  effect  on  the  antenna  impedance  near  these  resonances.  The  computed  reactance 
X  undergoes  seven  sign  changes  in  the  h/x  interval  between  0.03  and  0.1,  but  this 
could  not  be  shown  in  Fig.  6. 

The  antenna  impedance  near  the  resonance  peak  of  k^h  *  «  is  further  illus¬ 
trated  in  Fig.  7.  The  Impedance  is  computed  with  a  gradual  change  of  the  antenna 
length.  The  impedance  components  R  and  X  which  are  computed  based  on  an  assumed 

8  S 

sinusoidal  current  distribution  remain  nearly  constant,  but  the  impedance  component 
R  and  X  which  are  computed  using  a  two -term  trial  function  exhibit  significant 
changes.  The  two  sets  of  computations  give  nearly  the  same  impedance  figures  for 
Re  k0!  =  mt,  but  the  largest  differences  in  impedance  occur  near  Re  kg!  »  10. 5  n«, 
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where  the  antenna  reactance  is  changed  from  capacative  to  inductive.  However 
the  impedance  curves  are  not  symmetrical  with  respect  to  Re  k^i  *  10. 5  n*  and 
the  deviations  of  the  antenna  impedance  Z  from  Zg  values  are  decreased  in  magnitude 
with  an  increased  antenna  length,  lhe  amplitude  A  of  the  fundamental  sine  wave 
is  less  than  the  amplitude  B  of  the  surface  waves  in  the  vicinity  of  Re  k^i  - 
10. 5*.  However  the  surface  wave  have  considerably  less  effect  for  larger  values 
of  the  layer  thickness  h  and  for  longer  antenna  lengths  J,  when  Z  of  (35)  can 
be  expected  to  indicate  the  correct  order  of  magnitude. 


7.  SUMMARY  OF  CONCLUSIONS 

The  impedance  of  an  antenna  has  been  computed  in  the  presence  of  a  compressible 
plasma  by  considering  the  perturbations  of  the  antenna  current  distribution  by 
plasma  waves.  With  the  exception  of  very  thin  plasma  layers  the  antenna  Impedance 
can  be  determined  by  neglecting  the  presence  of  surface  waves.  The  impedance  of 
a  finite  antenna  was  shown  to  approach  the  impedance  in  the  presence  of  cold 
plasma  as  the  acoustic  velocity  is  decreased.  The  presence  of  an  Insulating 
layer  (or  of  an  ion  sheath)  around  the  antenna  also  tends  to  decrease  the  compressi¬ 
bility  effects. 

The  present  Impedance  computations  tend  to  justify  past  work  which  has  been 
carried  out  by  assuming  that  the  antenna  current  distribution  is  the  same  as  for 
cold  plasma.  On  the  other  hand,  calculations  which  neglect  the  presence  of  the 
electromagnetic  waves  and  attribute  the  antenna  current  variations  solely  to 
plasma  waves  [Cook  and  Edgar,  1966]  sure  obviously  in  error. 

This  behaviour  of  a  linear  antenna  is  quite  different  from  a  vide  slot 
antenna,  where  the  compressibility  of  the  plasma  effects  the  antenna  Impedance 
only  when  considering  the  surface  waves  in  a  direction  transverse  to  the  antenna 
aperture  [Galejs,  1966a].  The  waveguide  aperture  is  vide  in  terms  of  the  plasma 
wavelength  and  the  presence  of  a  compressible  plasma  would  have  negligible 


effects  on  the  impedance  of  a  linear  antenna  of  comparable  dimensions.  Further¬ 
more  the  surface  waves  in  the  linear  antenna  case  would  perturb  only  the  longi¬ 
tudinal  current  distribution  which  has  been  shova  to  have  small  effects  on  the 
antenna  impedance. 

The  antenna  of  the  present  investigations  was  idealized  as  a  current  sheath. 
For  thin  insulating  layers  the  surface  waves  supported  in  this  antenna  geometry 
are  different  from  Idle  waves  near  a  perfectly  conducting  antenna.  However  these 
differences  become  small  for  insulating  layers  of  a  thickness  of  the  order  of  the 
plasma  wavelengths  [Galejs,  1965a],  and  the  conclusions  of  the  present  analysis 
will  apply  also  to  metallic  antennas  if  the  ion  sheath  is  of  such  a  thickness. 
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Figure  2.  Antenna  Impedance.  Effects  of  Dielectric  Constant 
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Figure  3.  Antenna  Impedance.  Effects  of  Acoustic  Velocity. 
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Figure  6.  Antenna  Impedance  for  Plasma  Layers. 
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